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Abstract 



In this article we reconsider a version of quantum trajectory theory based on the stochastic 
Schrodinger equation with stochastic coefficients, which was mathematicaUy introduced in 
the '90s, and we develop it in order to describe the non Markovian evolution of a quantum sys- 
tem continuously measured and controlled thanks to a measurement based feedback. Indeed, 
realistic descriptions of a feedback loop have to include delay and thus need a non Markovian 
theory. The theory allows to put together non Markovian evolutions and measurements in 
continuous time in agreement with the modern axiomatic formulation of quantum mechanics. 
Oh To illustrate the possibilities of such a theory, we apply it to a two-level atom stimulated 

■Jii by a laser. We introduce closed loop control too, via the stimulating laser, with the aim to 

^ enhance the "squeezing" of the emitted light, or other typical quantum properties. Note that 

^ here we change the point of view with respect to the usual applications of control theory. 

In our model the "system" is the two-level atom, but we do not want to control its state, 
'— ' to bring the atom to a final target state. Our aim is to control the "Mandel Q-parameter" 

and the spectrum of the emitted light; in particular the spectrum is not a property at a 
^ single time, but involves a long interval of times (a Fourier transform of the autocorrelation 

function of the observed output is needed). 

00 

^ 1 Quantum trajectories and control 

Stochastic wave function methods for the description of open quantum systems are now 
widely used [l}|5] and are often referred to as quantum trajectory theory. These approaches are 
" . I very important for numerical simulations and allow the continuous measurement description 

J>" of detection schemes in quantum optics, namely direct, homodyne and heterodyne photo- 

detection (6]-[9]. In the Markovian case, the stochastic differential equations of the quantum 
trajectory theory can be interpreted in terms of measurements in continuous time because 



X 
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they can be related to positive operator valued measures and instruments 10 -12 , which are 
the objects representing observables and state changes in the modern axiomatic formulation 
of quantum mechanics. Moreover, these stochastic differential equations can be deduced 
from purely quantum evolution equations for the measured system coupled with a quantum 
environment, combined with a continuous monitoring of the environment itself. Such a 
representation is based on the use of quantum fields and quantum stochastic calculus |8||13|- 



18 



The whole quantum trajectory theory is well developed in the Markovian case, but to 
include memory effects is more and more important. A generalization by Diosi, Gisin and 
Strunz, based on the introduction of functional derivatives acting on the "past" inside the 
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stochastic Schrddinger equation (SSE), is able to describe dynamical memory effects 19 20 



but fails to have an interpretation in terms of continuous measurements 21 



Another way to include non Markovian effects and to maintain at the same time the 
continuous measurement interpretation is to start from the linear SSE and to generalize it 
by allowing for stochastic coefficients. This can be done without violating the axiomatic 
formulation of quantum mechanics and a non Markovian quantum trajectory theory can 



be developed in a mathematically consistent way 22 . The key point to get this result 
is that the non Markovian character is introduced by using stochastic coefficients in the 
linear version of the stochastic Schrodinger equation, without other modifications, and this 
allows to construct completely positive linear dynamics and instruments as required by the 
principles of a quantum theory. 

More recently, motivated by the growing interest for non Markovian evolutions, we begun 
to analyse possible physical applications of this theory. Some applications to systems affected 



by coloured noises and continuously monitored have been already developed 23 24 . This 
paper, instead, will be focused mainly on feedback control. 

In quantum optical systems, even when the Markov approximation for the reduced dy- 
namics is well justified, memory can enter into play when imperfections in the stimulating 
lasers are taken into account [25j and when feedback loops are introduced to control the 
system [4 26 - 33 . The so called closed loop control is based on the continuous monitoring of 
the system and, so, it fits well in the theory of measurements in continuous time. In some 
approximations, one can consider an instantaneous and very singular feedback and in this 
case the usual Markov framework is sufficient; however, more realistic descriptions of the 



feedback loop, including delay, need a non Markovian theory [9[|27 -33 



In this paper we present the non Markovian version of the theory of quantum meas- 
urements in continuous time, based on the SSE and the stochastic master equation (SME) 
with stochastic coefficients. Then, we develop such a theory, we explain how the stochastic 
coefficients can describe the most general measurement based feedback and also some imper- 
fections in the apparatus, and we show how to get the physical probabilities for the output 
of the observation, its moments and other related quantities of physical interest. Indeed, in 
quantum optical systems the moments of the stochastic output are connected to the Mandel 
Q-parameter and to the spectrum of the emitted light (homodyne and heterodyne spectra) 
and allow for the study of typical quantum properties of the emitted light, such as squeez- 



ing [9jj34j 35 . To illustrate these concepts we shall use a prototype model, a two-level atom 
stimulated by a laser, which is known to have a rich spectrum and to emit squeezed light 
under particular conditions. Our emphasis will be on the possibility of modelling a non per- 
fectly monochromatic and coherent stimulating laser and of modelling a measurement based 
feedback. We shall introduce closed loop control, via the stimulating laser, with the aim 
to control the squeezing in the observed spectrum, not to control the state of the system. 
Let us stress the change of point of view with respect to the usual applications of control 



theory 36 -39 . Here the "system" is the two-level atom, but we do not want to control 
its state, so as to bring the atom to a final target state, which is however possible inside 
our theory. Our aim is to control the properties of the emitted light; moreover, we want to 
control the spectrum, which is not a property at a single time, but involves a long interval 
of times (a Fourier transform in time is needed). 



2 The stochastic Schrodinger equation and the stochastic 
master equation 

The best way to introduce memory in quantum evolutions is to start from a dynamical 
equation in Hilbert space; this approach automatically guarantees the complete positivity of 
the evolution of the state (statistical operator) of the system. Moreover, considering the linear 
version of the SSE allows to construct the instruments related to the continuous monitoring 
even in the non Markov case [22]-p4]. We shall introduce first several mathematical objects. 
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from the linear SSE ([T]) to the instruments (14 1, and later, thanks to these latter, we shall 
give a consistent physical interpretation of the whole construction. 

Let !K be the Hilbert space of the quantum system of interest, a separable complex 
Hilbert space, with inner product (-I-). Moreover, we denote by the space of the 

bounded operators on by C the trace class and by §{3-C) C T(!K) the convex 

set of the statistical operators. 



2.1 The linear SSE and the reference probability 

Let us consider a reference probability space (il,3^, Q) with a filtration of cr-algebras {3^t)t 
satisfying the usual hypotheses, i.e. A g J with Q{A) = imphes A £ S'o, and = r\T>t -^t- 
As usual, we denote by uj the generic sample point in the sample space f2. The d + d' driving 
noises of the key SDE are defined in this filtered probability space; they are d continuous 
standard Wiener processes Bi, . . . ,Bd, and d' Poisson processes A^d+i, . . . , -/V^+rf'. Under the 
reference probability Q, all these processes are independent and are adapted, with increments 
independent from the past, with respect to the given filtration. Every Poisson process is 
taken with trajectories continuous from the right and with limits from the left (cadlag); let 
Aa: > be the intensity of Nk- 

We assume also we have a set of stochastic processes H{t), Li{t) (i = 1, . . . ,d + d'), with 
values in £(^K), such that H{t,uj)* = H{t,uj) and t Li{t), 1 1— H{t) are adapted processes 
with trajectories continuous from the left and with limits from the right (caglad, continuity 
in the strong operator topology). We assume also 

where M(T) does not depend on the sample point oj. 

Then, we introduce the linear SSE, for processes with values in 5f , 



d+d' 



d+d' 



i=l k=d+l 

d d+d' , J , s 



k=d+\ 



-lU(<_)d7Vfe(t), (1) 



with initial condition 



(/)(0) = e Jf, 



= 1. 



(2) 



Equation ([T]) is an Ito-type stochastic differential equation admitting a unique strong solution 
[22[ Proposition 2.1]. The solution is continuous from the right and with limits from the left 
(cadlag); the symbol t_ means to take the limit from the left. 

Here only bounded coefficients are considered, in order not to have mathematical com- 
plications, but generalizations to unbounded coefficients are of physical interest. Note that 
the filtration (J'j) can be taken bigger than the natural filtration generated by the driving 
noises B and A^, so that the processes L^, H can depend also on some other external noises. 



2.2 The linear stochastic master equation 

The SSE ^ can be translated into a stochastic differential equation for trace class operators. 
Indeed, adopting for a moment the Dirac notation, we can consider the rank one selfadjoint 
operator |0(t))((/)(t)| and compute its stochastic differential. In this way we get a closed hnear 
equation, which can be extended to a process cr(f) with values in T(!K) and whose initial 
condition can be any pure or mixed state. This is the linear stochastic master equation 



3 



(SME) (22) Propositions 3.2 and 3.4]: 

d 

da{t) = C{t)[(j{t^)]dt + ^(L,(i)a(t_) + a{t^)Ldty)dB,{t) 

J2 (^Lk{t)ait_)Lkitr -ait_)){dNkit)-\kdt), (3) 



1=1 

d' 



k=d+l 

with initial condition 



a{0)=poeHJi). (4) 
The operator £{t) is the stochastic LiouviUian 

d+d' d+d' 

C{t)[r] -i [H{t), r]~-Y^ {HtyUit), r} + J] U(t)TLdtr . (5) 

i=l i=l 

The SME ^ admits a unique strong solution. Typically, the solution (j(t) is not Markovian 
because C{t) depends on the past through the random operators H{t), Li{t). 

The propagator. In the following we shall need the fundamental solution, or propagator, 
A{t, s) of Eq. ([3]), i.e. the random linear map on T(IK) defined by a{s) cr{t). By construc- 
tion A{t, s) is completely positive and satisfies the composition law A{t, s) o A{s, r) = A{t, r) 
ioT < r < s < t. The propagator A{t, s) solves Eq. ^ with initial condition A{s, s) = 1. 

2.3 The new probability and the non linear equations 

Let us fix a non random state z G S(IK) and define the stochastic processes 

rr f ^^\\ 1+ \ \v[i,'^Y^o{i,^) ifp(t,w)7^0, 

P(0 — Tr{cr(i)}, p(t,:jj)-~\ u \ n W 

[z ifp(t,w) = 0, 

2ReTr{L,(t)p(i_)} , i = (7a) 
ik{t):=Tv{Lkit)*Lkit)p{t^)}, k = d+l,...,d + d'. (7b) 
By taking the trace of ([3|, we have that p{t) satisfies the Doleans equation 

d d+d' 
1=1 k=d+l ^ 

with p(0) = 1 and where the coefficients mi{t) and ifc(i) depend on the initial condition pQ 
2.3.1 The new probability 

The key property of quantum trajectory theory is that Eq. ([s]) implies that p{t) is a mean-one 
Q-martingale [22, Theorem 2.4, Section 3.1]. This allows us to define the new probabilities 



dp{t) - p(t_) <! J2 ^^mB^{t) + J2 ( ^ - 1 ) (dA^fe(i) - Afedt) \ (8) 



(F):^ [ p{T,Lj)QidLj)^EQ[p{T)lFl VF e Jt- (9) 
Jf 



Due to the martingale property of p{t), the probabilities Pj^ are consistent, in the sense that 
P*„(F) = P^JF) for F e J„ t > s > 0. 

Of course, the new probability Pj^ modifies the distributions of all the processes appearing 
in the theory, in particular the distribution of the processes Bi and Nk- 
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2.3.2 The Girsanov transformation 

A very important property is that a Girsanov- type theorem holds |22j Proposition 2.5, Re- 
marks 2.6 and 3.5]. Under Pj^, in the time interval [0,r], the processes 

Wj{t):^Bj{t)- f mj{s)ds, j^l,...,d, (10) 



^0 

turn out to be independent Wiener processes. Moreover, iVjj+i, . . . , iVjj+^' become simple reg- 
ular cadlag counting process of stochastic intensities id+i, ■ ■ ■ , id+d'', the meaning of stochastic 
intensity is given by the heuristic formula 

EpT^[dNk{t)\Jt] = ik{t)dt. 
From these results we have immediately 

EpT [B,{t)] = EpT [m,{t)] = 2ReEQ [Tr , 

I Ep. [Nk{t)] = EpT [zfc(i)] = Eq [Tr{Lfc(t)*Lfe(t)a(i_)}] . 
2.3.3 The non hnear SSE 

In this subsection let us consider again a pure initial condition po = \(t>o){4'o\ In this 
case we have a{t) — \(f>{t)) {(t>{t)\ and, so, 

p{t) = mt)\\\ m,{t) = 2Rcm_)\L,{t)^{t_)), i = 



where •ip{t) is defined by normalizing the random vector (j){t), solution of the linear SSE Q: 



m,uj)\\ if ||0(t,^)|| ^0, 

h if ||0(t,^)|| =0. 



Under the new probability Pj^ ([9]), the normalized Hilbert space process i{j{t) turns out 
to satisfy the non linear stochastic equation 22, Theorem 2.7] 



E f'^^^-^(i-))dA..W, (11) 



k=d+l 



with initial condition ■0(0) = 4'o ^-nd 

-. d-\-d' 1 ^ ^ ^ 1 f^+f^' 

kit):^-iH{t)~-Y,L,{trU{t)--Y,rn,{tf + -J2rn,{t)L,{t) + - ^ ^,{t) 

4=1 1=1 i=l k=d+l 



It is equation (111 which is usually called SSE and which is the starting point for powerful 



numerical simulations. 
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2.3.4 The non linear SME 



Going back to the case of a generic initial state Q , it is possible to show that the stochastic 
state p{t) defined by (|6| satisfies a non linear SME under the new probability Pj^ [22l Remark 
3.6]: 

d 

dp{t) = £{t)[p{t-)]dt + {U{t)p{t_) + p{t-)U{t)* - m,{t)p{t)) dW,{t) 

i=l 

+ E {^^Lk{t)p{t-)Lk{tr ~ p{t^)]{dNk{t)~ik{t)dt), (12) 

k=d+l ^^f'^^' ' 

/9(0) = po G S(^K); the operator Lit) is the stochastic Liouvillian ([5|. In [o] the reader can 
find a complete discussion of the relations among the four stochastic differential equations 
([T]), ([3]), ( [Tl] ), (12 1, in the purely diffusive case with deterministic coefficients. 



2.4 The continuous measurement process 

Let us introduce now the real processes 

d 



,,-1 "'0 



-ee,{t), e=l,...,mj, (13a) 



d+d' I, 

(0,t 



Ih{t) := / ahk{t~ s)dNk{s) +Wh{t), h ^ nij + 1, . . . ,mi + mj, (13b) 



where the integral kernels aej{t — s) and afik{t ~ s) are deterministic, and where ei(t) and 
Wh(t) can be stochastic processes. Even more general expressions than the ones in Eqs. 



( 13 1 could be considered. These processes will be interpreted as outputs of the continuous 



measurement, as explained in the next subsection. For this interpretation it is important 



the natural filtration generated by the processes (13), which we denote by (£4). We assume 



2.4.1 Instruments and a posteriori states 

Now, for t > 0, let us define the map- valued measure It [22] Remark 4.2]: 

It{E)[g]^ [ A{t,0,Lj)[g]Q{dio), e E-t, G T(Jf). (14) 

JE 

Such a measure has the properties: 

(i) VE G £t, 1-t{E) is a completely positive linear map on T(J{), 

(ii) \/g e T(J{), Va e £(J{), Ti {alt{-)[g]} is a-additive, 

(iii) Vg e T(J£), Tr{Itmg]} = Tr{rf. 

Such a map- valued measure It is called an instrument with value space {fl, £t) [l0 12 and 
it can be consistently interpreted as a quantum mechanical measurement on the system IK 
of the processes Jf(s) and Ih{s) in the time interval [0,t]; the instrument gives both the 
probability distribution of the output and the state changes conditional on the observation. 

According to the physical interpretation of the notion of instrument, the probability of 
the event i5 G £*, when the pre-measurement state is po, is given by 

Tr{It{E)[po]} - Eq [lETr{am - PUE) (15) 
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and this shows that the physical probabiHty for the observation of the output up to time t is 
indeed the one introduced in Eq. Q restricted to When t goes from to T, the family 
of instruments It gives a consistent description of a continuous measurement performed on 
the system. 



Then, Eq. ( 13 1 can be interpreted as the effect of the measuring apparatus which pro- 
cesses the ideal outputs Bi{t) and Nk{t) by the classical response functions aij and Uhk, and 
which degrades the outputs by adding some more noises ee{t) and Wh{t) due to the phys- 
ical realization of the apparatus itself. In quantum optics, the typical output current of an 



homodyne or heterodyne detector is of the form (13a I with 



a,j{t)^6,jFe{t), (16) 

where is the detector response function. In the final part of this section we shall see how 
to compute some relevant properties of the outputs under the physical probability. 

After all, we can say that (£t) j>o filtration generated by the outputs, i.e. £t contains 

only the events related to the observation of the outputs up to time t, while 3^t 1^ £-t, t > 0, 
is the filtration generated by all the processes involved in the continuous measurement and 
in the evolution, from the outputs to the unobservable noises. 

Now, let us take the conditional expectation on £t of the random state p{t) defined by 

^'--^'iJ'Wi^'i-iroii)^ (17) 

The interpretation is that p{t) is the conditional state one attributes to the system at time t 
having observed the trajectory of the output up to time t. Indeed p{t) is £t-measurable, thus 
depending only on the trajectories of the output in [0,<], and, moreover, one can directly 
check that [221 Remark 4.4] 

/ pit,L,)V'Jdoj)^ItiE)[po], yEeEf (18) 
Je 

The state p(t) is the a posteriori state at time t. 

In the extreme case £t — 3^4, which occurs for example when 3^t is generated by the 
processes Bi and Nk and just these processes are the observed output, we get p{t) = p{t). 
Therefore the evolution of the a posteriori state is completely defined by the non linear SME 
(12 1 satisfied by p{t), or, equivalently, by the SME ([3| for a(t), or even, as pure states are 
mapped to pure states, by the SSE ([T]) for (j){t) or the non linear SSE (11) for ip{t). Let us 
stress that the evolution is not Markovian due to the randomness of the operator coefficients 
L,{t), H{t). 

When £( C j'j^ the a posteriori state p{t) has a non Markovian evolution which typically 
does not even satisfy a differential equation. In this case the SSE ^ and the SME ^ have 
to be interpreted as an ideal unravelling of the physical evolution of p{t) which allows us to 
consistently define it, by ([3| and ([Tt]), and allows us to compute, at least numerically, all the 
quantities of physical interest (that is to define the instruments It). 

The randomness of the coefficients Li{t) and H{t) and the distinction between £t and 
3^t provide the SME ^ with a great flexibility which allows us to model non Markovian 
features of the dynamics due both to some environmental noises and to measurement based 
feedback loops. 

2.4.2 A priori states 

When the output of the continuous measurement is not taken into account, the state of the 
system at time t is given by the mean state 

r,{t) :=Ep. [p(i)] =EQ[a(t)] =EQ[^(t,0)[po]] =It(r!)[po]. 
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The state ri{t) is the a priori state at time t. Note that EQ[yt(t, 0)[-]] is a completely pos- 
itive, trace preserving, linear map, i.e. a quantum channel in the terminology of quantum 
information. 

From the SME ^ we get 

||^77(t) =Eq [C{t)[a{t)]] = I^C{t,u;)[a{t,u;)MdLu), (19) 

which is not a closed differential equation when C{t) is stochastic, contrary to the Markov case 
[9| Section 3.4]. By the projection operator technique a closed integro-differential equation 
for the a priori state ri{t) could be obtained |24] (an evolution equation with memory), but 
this equation is too involved to be of practical use. Again 

and ^ are an unravelling of 
a non Markovian evolution. By the non linear SME ([T2| we have also 



dt 



0<t<T. 



2.4.3 Spectra and moments of the diffusive outputs 



Let us consider an output current Je as given by Eq. (13a I; Jg is a diffusive stochastic 
process and its spectrum is given by the classical notion [40|. If Ji is at least asymptotically 
stationary and the limit in ( 20 ) exists in the sense of distributions in /i, then the spectrum 
of Je is defined by 

i-T 



lim — EpT 



Mt)dt 



(20) 



In the case (16 1 without extra noise, ee{t) = 0, and for a detector response function going to 
a Dirac delta, that is in the limit case Je — Bg, the spectrum becomes 



lim — EpT 

T-^+oo T "0 



e'^*dS^(i) 



(21) 



The spectrum depends on the distribution of the current Je{t), which is the output of 
the continuous measurement on the system, and, so, Se gives information on the monitored 
system. Actually, the spectrum gives information also on the carrier which mediates the 
measurement. For example, in quantum optics, the system is typically a photoemissive source 
which is continuously monitored by homodyne or heterodyne detection of its emitted light. 
In this case, Se gives information both on the photoemissive source and on the fluorescence 
light. For example, Se < I reveals squeezing of the emitted light. 

An expression for the autocorrelation function needed in the computation of the spectrum 
can be obtained by generalizing the techniques used in the Markovian case [9] Section 4.5]. 
When the Liouville operator ^ and Be are independent (which implies that Be is not used 
for the feedback), we get 



dtds 



EpT [Be{t)Be{s)]=5{t- s) + be(t,s)+be{s,t), 



(22a) 



he{t,s) = l(o,+oo)(i - s) Ej; 



Tr<i {Le{t) + Le{t)*) 

X A{t-,s) [Le{s)(j{s^) + (7{s^)Le{sy 



(22b) 



Let us set ne{t) = Ep* [m£(t)] and assume that the limit Uoo '■— lini(->.+oo ^^^(i) exists. 
Then, we obtain the decomposition of the spectrum in the elastic part and the inelastic one 
(the spectrum of the fluctuations) 

Se (m) = '5'gI (m) + '5'inci ip) , Sciin) = 27r (5 (p) , (23a) 
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SinoKA*) = 1 + lim I; / dt [ ds COS iJ,{t- s)de{t,s), (23b) 

de{t,s) := bt{t,s) - l(o,+oo)(^- s)nryoni{t) = l(o,+oo)(^ - «) 

X Eq [Tr {(L,(i) + Li{ty)A{t^,s) [L,(s)a(.s_) + a(s_)L,(s)* - n^a{s^)]}] . 

2.4.4 The Mandel Q-parameter of the counting outputs 

When we consider direct detection in quantum optics, in the ideal case of noiseless counter, 
the output of the measurement is one of the counting processes, say Ik = N^, which means 



to take ahk(t) = ^hk and Wh(t) = in Eq. (13b I. In this case a typical quantity is the Mandel 
Q-parameter, defined by 

Ep/[7V,(io + <)-iVfc(to)] ^' 

As for a Poisson process this parameter is zero, in quantum optics it is usual to say that in 
the case of a positive Q parameter one has super-Poissonian light and sub-Poissonian light 
in the other case. Sub-Poissonian light is considered an indication of non-classical effects. 

In quantum trajectory theory one can find expressions for the moments also in the count- 
ing case and we get 

Q.(i;^o) = |^|^-A4(i;M, 

Mk{t\ to) 

/•to+t 

Mk{t;to) :=EpT [Nk{to + t) ~ Nk{to)] = / dsEQ [Tr {Lfc(s)*Lfc(s)CT(s_)}] , 

™ Jto 

Vkit;to) :=EpT ]^{Nk{to+t) - 7Vfe(<o))'] - Mk{t]h). 

When Lk{s)* Lk{s)A{s-^r) is Q-independent from Nk{r), which happens when Nk is not 
used for feedback, we get 

Vkit;to)^2 ds f drEQ[Tr{Lfc(s)*Lfc(s)^(s_,r)[Lfc(r)a(r_)L,(r)*]}]. 

Jto J to 

3 A paradigmatic model: the two-level atom 

As an application of the theory we consider a two-level atom stimulated by a laser; it is an 
ideal example, but is sufficiently rich and flexible to illustrate the possibilities of the theory. 
The Hilbert space is and the Hamiltonian part of the dynamics is given by 

H{t)=Ho + Hf{t), Ho^'ja,, H fit) =m + f{t)a+, (24) 

where t'o > is the resonance frequency of the atom and CTz, a± are the usual Pauli matrices. 
The function / is the laser wave, which can be noisy and can be controlled by the experi- 
menter; a concrete choice for / is given below. 

Let us complete the model by choosing the noise-driven terms in the SME (Is]), which we 
call channels in the sequel. A sketch of the ideal apparatus is given in Figure 

We consider two diffusive channels realized by heterodyne or homodyne detectors of the 
emitted light with local oscillators represented by the functions hj : 



d = 2, Lj{t) = hj{t)aja_, \hj{t)\ = 1, G C. 

The observation of the light in channel 1 will be used to control the stimulating laser light. 
The light in channel 2 is only observed, in order to analyse its properties. We are interested 
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^3(t) 


photo- 




-counter 



lost light 

) channel 4 



sp. analyser 










Figure 1: The electro- modulator EOM acts as a phase modulator. The beam-splitter near EOM 
has transmissivity and reflectivity nearly equal. The second beam-splitter has a transmissivity 
much smaller than the reflectivity. In this way the two local oscillators are much more intense 
than the light stimulating the atom, as required by the homodyne configuration. 



in controlling the properties of this part of the emitted light by controlling the atom via the 
stimulating laser; in particular we are interested in generating squeezed light. 

In the case of homodyne detection the local oscillator is fed by the light produced by the 
stimulating laser and we have hj{t) — e~'^J /fft)\ • This is the case illustrated in Figure 1 

In the case of heterodyne detection the local oscillator is fed by an independent laser. For 
this light we could use for instance the so called phase diffusion model, 

hj{t) — exp {— icj — Wjt — i/c_ji?_j(i)} , 

where B„i and i?_2 are extra noises (independent standard Wiener processes). The function 
hj represents a nearly monochromatic wave with a Lorentzian spectrum centred on Vj . 
We introduce also four jump channels: 

2 4 

d'-4, Lk{t)=Lk, 13k eC, 7>0, n>0, ^ |a,|' + ^ = 7, 

i=l k=3 

Lk = Pk(y-, fc = 3,4, 5 = /36cr+, \P5? = = in. 

The jump channels 3 and 4 are electromagnetic channels: channel 3 represents the emitted 
light reaching a photo-counter {direct detection)^ while channel 4 represents the lost light, 
which is not observed. The output of channel 3 could be used again as a possible signal for 
closed loop control, but we use it here only to see properties of direct detection. As channel 
2 is used to detect properties of the field quadratures (squeezing), the analogous channel 3 
is used to explore the counting statistics of the emitted photons. 

The counting channels 5 and 6 are used only to introduce dissipation due to a thermal 
bath; these channels are not connected to observation. So, channels 4, 5, 6 are introduced 
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only to give dissipative effects in the system dynamics; to represent them in the jump form 
(as done here) or in the diffusive form has no influence on the physical quantities. 

The stimulating laser light can be noisy and can be controlled by the output of the 
diffusive channel 1. In mathematical terms, / can be an adapted functional of Bq (extra- 
noise) and Bi (feedback of the output of the homodyne detector 1). 

For this model we need a probability space (fi,?, Q) where _B_2,...,i?2 are Wiener 
processes, N^, . . . , Nq are Poisson processes and they are all independent. According to the 
notations of Section [2j the filtration (9^() is generated by all these processes, while (£j) is 
generated by Bi, B2 and N^. 

From Eq. ([s]) and the assumptions of the model we get the random Liouville operator 

C{t)[T] = — i [Hq + Hf{t), t] + jcr^ra^ — - {cr+a^^T} + ■yn {u^ra^ + Oj^to^ — t) . 

Note that the whole of the randomness is in the wave f{t) and it is due purely to noise in 
the laser light and to feedback, but this is enough to have a non Markovian evolution of the 
atom. 



3.1 The laser wave 

Let us denote by fa{t) the laser wave without feedback. We assume the laser to be not 
perfectly coherent and monochromatic and we describe this again by the phase diffusion 
model, i.e. 



-i[i?+i't+fcoBo(t)] 



i?G(-7r,^]; (25) 



rifi is known as Rabi frequency and 
spectrum of this wave; we get 



— v as detuning. It is easy to compute the 



hm - E 



JfJ.t 



fomt 



O 21.2 



fc4 + 4(^-i.)^ 



So, /o represents a wave of constant amplitude |/o(^)l^ — ^r/^ ^-nd Lorentzian spectrum. 
The effect of such a stimulating laser wave in the heterodyne spectrum of a two-level atom 
has been studied in Refs. 



18 25 



In that papers, however, the presentation of the theory 
of continuous measurements was based on the use of Bose fields and quantum stochastic 
calculus. 

The present theory gives also the possibility of describing feedback control, by assuming 
that f{t) is a functional of the output current Ji(s), < s < t. We take Ji to have the form 
given in Eqs. (13a), (16 1 with a very simple detector response function; to be concrete we 
take 



Ji(t) = k 



fits) 



dBi(s), 



k e 



>c> 0. 



(26) 



The functional dependence of / on Ji is determined by the implementation of the feedback 
mechanism, which should be chosen according to the aims of the controller and to the physical 
feasibility. Then, in principle, the quantities of interest could be computed by starting from 
simulations of the non linear SSE of the model. However, in order to have an analytically 
computable spectrum and to have an idea of the possible behaviours of the model, we consider 
very simple forms for /. 



A first choice is to modify the wave ( 25 ) by a simple linear amplitude modulation and to 
take 



f{t) ^ ie-'[^+'^*+'="^°(*)i + 



ce 



ii5 



'Mt-6)] 



c e 



e (- 



(27) 
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S > represents a possible delay. By taking S = 0, k = ^ x/{2'it), k t +00 one formally 
obtains J\ (t) = Bi (t) , the limiting case of the "Markovian feedback" of Wiseman and Mil- 
burn [27{[28] ; in [35] and in [9| Chapt. 10] this feedback scheme has been formulated in terms 
of SMEs and applied to the homodyne spectrum of a two-level atom. The effect of delay in 
some models has been studied in 28j 33 4l] . 

A second choice of the control is to introduce a phase modulation of the wave ( 25 ) . By 
taking this phase modulation linear in the output current we get 



fit) 



R ^~i[i)+ut+koBoit)+kiJi(t-S)] 



S>0, he 



(28) 



Our aim is only to show that it is possible to construct in a consistent way a theory with a 
non Markovian feedback and to demonstrate that, even in simple cases, the feedback can be 
used to control the squeezing. So, we simplify further on the expression ([28| and we take 
no delay, S — 0, and x — 0, k = 1. Then, we have Ji{t) = Bi{t); note that the formal 
instantaneous output is Bi{t) and, thus, Bi{t) is the integrated current: f{t) depends on 
the whole past of the signal detected by the homodyne detector 1. Our final choice for the 
laser wave is 



fit) 



u{t) = d + ut + koBo{t) + kiBi{t). 



(29) 



Whatever form we take for /, the feedback acts on the stimulating light and on the local 
oscillators, as illustrated by Figure [T] As / enters the Hamiltonian (24 1, we can control to 
some extent the dynamics of the atom. This possibility can be used to control the system 
state or the time to reach a final state [30{[4l] . Instead, what we want to study in this paper 
is how the feedback can be used to control the output of the quantum system, precisely the 
squeezing in channel 2 or the Mandel Q-parameter in channel 3. 



3.2 Control of the homodyne spectrum 



According to Eqs. ( 23 1, to compute the spectrum of the light in channel 2 we need to compute 



first the quantities 71.2 (i) £^nd d2it,s). The best way is to introduce a kind of stochastic 
rotating frame, that is to make a unitary transformation and to define the quantities 



A(t,s)[r] :=e5"(*)^M(t,s) 



e 2 



u{t)(7z 



where A{t, s) is the propagator defined at the end of Section 2|2.2 and 



(30) 



By using ([s]), ([29| and computing the stochastic differential of ^[t) we obtain 

2 

de(t) -£(t)K(L)]di + ^A(t)K(t-)]di?.(t) 



i=0 



(^Lk4it-)Lk* -iit-)) {dNkit)-\kdt), (31) 



fc=3 



* 1 ^ 

^t)[T] = - - [Ai^a^ +nji(T^, t] +7(t_tct+ - -{-P+,r} 



ifci (9iit)P+T(T+ - gi{t) a-TP+^ + jn {a-Ta+ + a+ra^ - r) , (32) 
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^0 M = ^ ^0 [ctz , t] , Vi (t) [t] = gi (t) a_T + gi {t)T(j+ + ^ [^^ , r] , 



+ 



By the hnk ( [30| ) between ^ and A, we have that also A(t, s) satisfies Eq. (31) with initial 
condition A(s, s) = 1. 
Then, we get 

n2{t) = EpT [m2(t)] = 2ReEQ [Tr {L2(t)fT(L)}] = Eq [Tr{2?2(i) , 

d2{t,s) = l(o,+oo)(i - s)Eq [Tr{2?2(t) [A(L,s) [2?2(s) - ^00^(5-)]]}] . 

All terms inside the expectation values are random and in general 712 and d2 can not be 
computed explicitly. 

3.2.1 Homodyning in channels 1 and 2 

In the case of homodyning as in Figure [ij we have Qjit) — |aj|e^"'^, j — 1,2, which gives 



I?2, non random and independent of time; the operators in (31 ) were already non 
random. This is the key point which makes the model computable. By Eq. ( [sT] ), A(t_, s) and 
^(s_) turn out to be Q-independent and we get EQ[A(t_, s)] = e^(*^*\ Eq[^(s_)] = e^** [po]- 
Finally, we get 

n2{t) = 2|a2|ReTr(e"'''V+e^* [pA , = lim n2{t), 



d2{t,s)= 2|a2|^ReTr|e~' 



e'^V_e^^ [po]+h.c.]} 



n^n2{t) 1(0, +00) - s), 



where h.c. means "hermitian conjugate term". The modified Liouville operator C can be 
written in the Lindblad form as 



C[r] 



H. T 



1 



+ (7-|«il') U-r<J+--{P+,T}]+KTK* 



1 k 

- - {K*K, '^}+ ^ {a-zTcr^ ~ t) + jn {(T-T(t+ + a+ra- - r) , (33) 



|ai|e'^i(7- 



By Bloch equation techniques we can compute the homodyne spectrum ( 23 ) of the light in 
channel 2 [oj Part II] . One can show that the Liouville operator ( 33 ) has a unique equilibrium 
state and that, because of the limits contained in (23), any dependence on the initial state po 
disappears from the expression of the spectrum. The final result is 5'2(/i) = Sc\{p) + S'inoi(A*): 



Sc\{p) = 27r |q!2|^ v^SifJ'), V cos?92 di + sin 1^2 ^2, 



'5'inci(Ai) = 1 + (cos 1^2 sm'&2 0) 



A = 



-Av ^ fiij; + fci |ai I cos 7?i 







2 




COS 1^2(1 + d3)\ 
A ( sin 792(1 + ^3) +VU 

^— fci jail sin??!^ 

U = ( fci |q!i| COS"!?! 

7 



(34) 
(35) 

(36) 
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d = -A- 



u, 



r = (2n + l)7- 



ki- 



rn 



Now we can study the dependence of this spectrum on our parameters, feedback included, 
and learn how to control the emitted light. By certain choices of the parameters one can 
obtain Siacl{^J■) < Ij which is interpreted as squeezing of the light in the channel 2. 

As an example, let us take fco = 0, 7 = 1, n = 0, |ai|^ = \a2f = 0.45. We choose 
the other parameters to have the deepest well for S'inoi in a given point /i. We denote by 
Sol = 27r |q!2|^ the coefficient of the delta-spike in the elastic part. 



-1.5708 we get Sei = 1.0245 and 
-0.1540 we get 



• For = 0.366, Av = 0, fci = 0.3371, = -tt, 1S2 
a minimum in /i = of S'inoi (0) = 0.8172. 

• For Qr = 1.6150, Aj/ = 1.3833, ki = 0.3213, = -1.9307, ^2 
Sol = 1.4214 and a minimum in fi ~ 2 of S'inoi(2) = 0.8621. 

• For ^Ir = 3.1708, = 2.5576, ki = 0.3249, i?i = -1.7863, i?2 = -0.0760 we get 
Sol — 1.5356 and a minimum in fi — A of S'inci(4) — 0.8572. 

In Figure [2] we plot the inelastic spectrum as a function oi x — fi for a case with control 
and a case without control. The parameters are optimized to have the deepest minimum in 
/X = 2. 




Figure 2 

I i2 I i2 
\Oil\ = 02 



n 



R 



Squeezing control. S'inoi (/^) with and without feedback for 7 = 1, /cq = = 0, 
0.45 and: (solid line) ki = 0.3213, ??i = -1.9307, ^2 = -0.1540, Au = 1.3833, 
1.6150; (dotted line) ki = 0, t?2 = -0.1784, Au = 1.4937, Qr = 1.4360. 



The values of the minima reported above are very similar to the ones found in [9, p. 244], 
where a very different scheme of feedback was considered (the instantaneous and singular 
feedback a la Wiseman-Milburn 27 discussed after Eq. (|27])). The two schemes are very 



different: one is a limiting case of amplitude modulation and the other an over-simplified 
case of phase modulation; it is not obvious a priori that one can obtain similar values for 
the optimal minima in the two cases. The difference of the two schemes is stressed also 
by the different values of the other control parameters which realize the minima in the two 
cases. For instance, the best minimum S'inoi (4) = 0.8572 in /i = 4 in the present case is 
obtained for the Rabi frequency ^Ir = 3.1708 and the detuning Ah' — 2.5576, while in the 
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case of [9j Sect. 10.5.3] the best minimum S'inei(4) — 0.8565 is obtained for fifl = 3.0721 and 
Az/ = 2.6310. Without any feedback the best minimum is S'inci(4) = 0.8830 for = 2.9001 
and Ai/ = 2.8077. 

By using quantum fields and quantum stochastic calculus it is possible to prove that the 
Heisenberg uncertainty principle implies that the product of S'inoi(M)j computed for a certain 



value times the same quantity for t?2 + 7i'/2 is always not less than 1 34 



^inollM; '&2) S'inol(Ai; 1?2 + 7r/2) > 1. 

In Figure [3] we plot the previous case with feedback, taken with three different choices of 1^2 
in order to see the dependence on -di, which changes the field quadrature under monitoring. 
We plot a case with a pronounced squeezing, the case with a 7r/2 shift in (anti-squeezing) 
and the case with a shift of 7r/4, just to see the behaviour of an intermediate field quadrature. 



A 




Figure 3: Effect of Heisenberg uncertainty on S'incK/t^)- The parameters are 7 = 1, A;o = n = 0, 
\axf = \a2\^ = 0.45, ki = 0.3213, t?i = -1.9307, Ai/ = 1.3833, = 1.6150 and: (solid line) 



-0.1540; (dotted line) ^2 



0.1540; (dashed line) t?2 



0.1540. 



3.3 Control of the photon statistics 

We maintain the same feedback as before, but now we study the Mandel Q-parameter of the 
direct detection = N^. We consider only the stationary regime at long times and we take 



Qz{t) 



lim Qz{t;h) 



hm ("iM^M^MalMo; 



By using again Bloch equation techniques to compute the expressions of Section [2|2. 4. 4[ we 
get 

hnr M3(i;to) = ^ l/Sgl' (1 + dg) , 

to— i-+oo Z 



Qsit) = lAI 
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where d, A are defined in Eqs. (36), (37). By taking a long time interval we get 



Q3:= lim QsW^-l/^sl' 

r— f +00 

According to the choice of the parameters we can obtain a sub-Poissonian or a super- 
Poissonian Q-parameter: fco = 0, 7 = 1, n = 0, |ai|^ = |/33|^ = 0.45, 

• Ai/ = 0, /ci = 1.0126, nji. = 1.0063, = tt: Qg = -0.5094; 

• Av^O, ki=0,QR = 0.7071: Q3 = -0.3375; 

• Av^2,ki^ 2.8515, = 2.3516, z9i = 2.6914: Q3 = -0.4356; 

• Ai/ = 2, /ci = 0, Qr = 2.9155: Q3 = 0.0860. 

Moreover, one can see numerically that in the case of no feedback, fci = 0, and Ai/ = 2, we 
have Q3 > for all ^Ir > 0. So, in this case feedback control is essential to get sub-Poissonian 
light. 

Finally, let us observe that channels 2 and 3 are completely analogous; only the two 
detection apparatuses are different. What comes out from the numerical computations is 
that no relation seems to hold between Qa < and squeezing; with the parameters which give 
squeezing we have Q3 = 0.0602 (for ^r = 1.6150, fci = 0.3213, i?i = -1.9307, Av = 1.3833) 
and Q3 = 0.09508 (for 17^ = 3.1708, ki = 0.3249, i9i = -1.7863, Av = 2.5576). 




4 Conclusions 

In this article we have shown how the SSE or the SME with stochastic coefficients allows 
to model a wide class of non Markovian quantum evolutions in agreement with the axioms 
of quantum mechanics. Let us stress that the theory is fully "quantum" in that it gives 
completely positive evolutions and instruments and it is able to describe quantum effects as 
squeezing and anti-squeezing. This statement can be strengthened by showing, at least in the 



Markov case 17 18 , the equivalence with a theory containing quantum systems, quantum 
fields, unitary evolutions and observation of compatible field observables. 

What we have shown in the non Markov case is that our approach can describe several 
evolutions with memory effects, both in the case of open but unobserved systems and in the 
case of continuously measured systems. 

The fundamental ingredient is to start from the linear SSE ([T]) : even if the coefficients are 
stochastic, starting from a linear equation for pure (unnormalized) states ensures a linear, 
trace preserving and completely positive mean evolution, which therefore can describe the 
evolution of an open system. Typically, such evolution is not Markovian, as the coefficients 
are stochastic, and therefore it does not satisfy a master equation, but it is defined by a SSE 
(or, equivalently, by a SME) plus an expectation. 

Moreover, the unravelling of this evolution provided by the SSE and the SME can have a 
physical interpretation in the theory of continuous measurements. Besides the linearity of the 
SSE ([!]) and the stochasticity of its coefficients, the other key ingredient is the introduction 



of the outputs processes Je and Ih ( 13 ) and of their natural filtration £t eventually smaller 
than the initial filtration 3^t' combining a linear equation for pure (unnormalized) states, 
stochastic coefficients and conditional expectations gives instruments It which can model a 
non Markovian continuous measurement. Indeed, it is possible to describe not only several 
sources of classical noise, such as imperfections of the apparatus or coloured noises, but also 
the most general measurement based feedback. Thus we can handle also feedback schemes 
with some delay in the feedback loop. Again, both the evolutions of the a posteriori and 
a priori states do not satisfy directly closed equations, but they are defined by stochastic 
differential equations plus (conditional) expectations. 

Moreover, as we are considering continuous measurements, our interest is not only for the 
evolution of the observed system, but also for the outputs of the observation: the outputs 
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are quite general, they have both diffusive and counting eomponcnts, and we study their 
distribution, providing useful formulae for the moments and the spectra of the diffusive 
outputs and for the Mandel Q-parameters of the counting outputs. 

To present an application of the theory, we consider a two-level atom, stimulated by a non 
perfectly monochromatic and coherent laser, and we introduce a thermal bath, a continuous 
measurement by detecting its fluorescence light, and a closed loop control by modulating the 
stimulating laser according to the observed output. The resulting evolution of the atom is 
not Markovian because of the imperfections of the laser and because of the chosen feedback 
loop. 

More precisely, the fluorescence light is divided in several channels where detectors are 
positioned; one of these performs homodyne detection and its output is used in the feedback 
loop with the aim of controlling, not the atom itself, but the light emitted in the other 
channels. The chosen feedback provides a non Markovian example where calculations can 
be explicitly performed up to find the homodyne spectrum and the Mandel Q-parameter of 
the fluorescence light and their dependence on the control parameters. Thus we study the 
role of the feedback loop in controlling the qiiantum properties of the emitted light, such 
as squeezing and sub-Poissonian statistics, we show that these quantum properties can be 
enhanced by feedback, but we find that there is no relation between them. 

Even keeping our interest confined to the two-level atom, several further developments 
are possible, such as using the feedback here introduced to control properties of the atom 
itself or other properties of the emitted light, maybe local in time, or such as introducing 
other measiirement and feedback schemes and then studying (numerically) the optimal one 
associated to a given target (e.g. "maximal squeezing" of the fluorescence light). 

References 

[1] Carmichael, H. J. 1993 An Open System Approach to Quantum Optics. Lectures Notes 
in Physics m 18. Berlin: Springer. 

[2] Zoller, P. & Gardiner, C. W. 1997 Quantum noise in quantum optics: the stochastic 

Schrodinger equation. In Fluctuations quantiques. (Les Houches 1995) (ed. S. Reynaud, 
E. Giacobino & J. Zinn- Justin), pp. 79-136. Amsterdam: North-Holland. 

[3] Breuer, H.-P. & Petruccione, F. 2002 The Theory of Open Quantum Systems. Oxford: 
Oxford University Press. 

[4] Gardiner, C. W. & Zoller, P. 2004 Quantum Noise: A Handbook of Markovian and Non- 
Markovian Quantum Stochastic Methods with Applications to Quantum Optics, Springer 
Series in Synergetics. Berlin: Springer. 

[5] Carmichael, H. J. 2008 Statistical Methods in Quantum Optics, Vol 2. Berlin: Springer. 

[6] Barchielli, A. 1990 Direct and heterodyne detection and other applications of quantum 
stochastic calculus to quantum optics. Quantum Opt. 2, 423-441. 

[7] Barchielli, A. & Belavkin, V. P. 1991 Measurements continuous in time and a posteriori 
states in quantum mechanics. J. Phys. A: Math. Gen. 24, 1495-1514. 

[8] Holevo, A. S. 2001 Statistical Structure of Quantum Theory Springer Lecture Notes in 
Physics m 67. Berlin: Springer. 

[9] Barchielli, A. & Gregoratti, M. 2009 Quantum trajectories and measurements in con- 
tinuous time: the diffusive case. Lecture Notes in Physics 782. Berlin: Springer. 

[10] Davies, E. B. & Lewis, J. T. 1970 An operational approach to quantum probability. 
Commun. Math. Phys. 17, 239-260. 

[11] Davies, E. B. 1976 Quantum Theory of Open Systems. London: Academic Press. 

[12] Ozawa, M. 1984 Quantum measuring processes of continuous observables, J. Math. 
Phys. 25, 79-87. 



17 



[13] Barchiclli, A. & Lupicri, G. 1985 Quantum stochastic calculus, operation valued 
stochastic processes and continual measurements in quantum mechanics. J. Math. Phys. 
26, 2222-2230. 

[14] Barchielli, A. 1986 Measurement theory and stochastic differential equations in quantum 
mechanics. Phys. Rev. A 34, 1642-1649. 

[15] Belavkin, V. P. 1988 Nondemolition measurements, nonlinear filtering and dynamic 
programming of quantum stochastic processes. In Modelling and Control of Systems 
(ed. A. Blaquiere). Lecture Notes in Control and Information Sciences 121, pp. 245- 
265. Berlin: Springer. 

[16] Belavkin, V. P. 1989 A new wave equation for a continuous nondemolition measurement. 
Phys. Lett. A 140, 355-358. 

[17] Barchielli, A. & Paganoni, A. M. 1996 Detection theory in quantum optics: stochastic 
representation. Quantum Semiclass. Opt. 8, 133-156. 

[18] Barchielli, A. 2006 Continual measurements in quantum mechanics and quantum 
stochastic calculus. In Open quantum systems III (ed. S. Attal, A. Joye & C.-A. Pillet). 
Lecture Notes in Mathematics 1882, pp. 207-291. Berlin: Springer. 

[19] Diosi, L. & Strunz, W.T. 1997 The non-Markovian stochastic Schrodinger equation for 
open systems. Phys. Lett. A 235, 569-573. 

[20] Diosi, L., Gisin, N. & Strunz, W.T. 1998 Non-Markovian quantum state diffusion. Phys. 
Rev. A 58, 1699-1712. 

[21] Wiseman H. M. & Gambetta J. M. 2008 Pure-state quantum trajectories for general 
non-Markovian systems do not exist. Phys. Rev. Lett. 101, 140401. 

[22] Barchielli, A. & Holevo, A. S. 1995 Constructing quantum measurement processes via 
classical stochastic calculus. Stoch. Proc. Appl. 58, 293-317. 

[23] Barchielli, A., Pellegrini, C. & Petruccione, F. 2010 Stochastic Schrodinger equations 
with coloured noise. EPL 91, 24001. 

[24] BarchieUi, A., Di Telia, P., Pellegrini, C. & Petruccione, F. 2011 Stochastic Schrodinger 
equations and memory. In Quantum probability and related topics (ed. R. Rebolledo & 
M. Orszag), QP-PQ: Quantum Probability and White Noise Analysis 27, pp. 52-67. 

Singapore: World Scientific. 

[25] Barchielli, A. & Pero, N. 2002 A quantum stochastic approach to the spectrum of a 
two-level atom. J. Opt. B: Quantum Semiclass. Opt. 4, 272-282. 

[26] Belavkin, V. P. 1983 Theory of the control of observable quantum systems. Automat. 

Remote Control 44, 178-188. 

[27] Wiseman, H. M. & Milburn, G. J. 1993 Quantum theory of optical feedback via homo- 
dyne detection. Phys. Rev. Lett. 70, 548-551. 

[28] Wiseman, H. M. & Milburn, G. J. 1994 Squeezing via feedback. Phys. Rev. A 49, 
1350-1366. 

[29] Giovannetti, V., Tombcsi, P. & Vitali, D. 1999 Non-Markovian quantum feedback from 
homodyne measurements: the effect of a non-zero feedback delay time. Phys. Rev. A 
60, 1549-1561. 

[30] Wang, J. & Wiseman, H. M. 2001 Feedback-stabilization of an arbitrary pure state of 
a two-level atom. Phys. Rev. A 64, 063810. 

[31] Wang, J., Wiseman, H. M. & Milburn, G. J. 2001 Non-Markovian homodyne-mediated 
feedback on a two-level atom: a quantum trajectory treatment. Chem. Phys. 268, 221- 
235. 

[32] Nishio, K., Kashima, K. & Imura, J. 2009 Effects of time delay in feedback control of 
linear quantum systems. Phys. Rev. A 79, 062105. 



18 



[33] Wiseman H. M. & Milburn G. J. 2010 Quantum Measurement and Control. Cambridge: 

Cambridge University Press. 

[34] Barchielli, A. & Grcgoratti, M. 2008 Quantum continual measurements: the spectrum 
of the output. In Quantum probability and related topics (ed. J. C. Garcia, R. Quezada 
& S. B. Sontz). QP-PQ: Quantum Probability and White Noise Analysis, 23, pp. 63-76. 
Singapore: World Scientific. 

[35] Barchielli, A., Gregoratti, M. & Licciardo, M. 2009 Feedback control of the fluorescence 
hght squeezing. EPL 85, 14006. 

[36] Belavkin, V. P. & Edwards, E. 2008 Quantum filtering and optimal control. In Quantum 
Stochastic and Information (ed. V. P. Belavkin & M. Gu^a) pp. 143-205. Singapore: 
World Scientific. 

[37] Boutcn, L. & Van Handel, R. 2008 On the separation principle in quantum control. 
In Quantum Stochastic and Information (ed. V. P. Belavkin & M. Gu\,a) pp. 206-238. 
Singapore: World Scientific. 

[38] Gough, J. 2008 Optimal quantum feedback for canonical observables. In Quantum 
Stochastic and Information (ed. V. P. Belavkin & M. Gu^a) pp. 262-279. Singapore: 

World Scientific. 

[39] James, M. R. 2008 Feedback control of quantum systems. In Quantum Stochastic and 
Information (ed. V. P. Belavkin & M. Gu^a) pp. 280-299. Singapore: World Scientific. 

[40] Howard, R. M. 2002 Principles of random signal analysis and low noise design, the 

power spectral density and its applications. New York: Wiley. 

[41] Combes, J. & Wiseman, H. M. 2011 Quantum feedback for rapid state preparation in 
the presence of control imperfections. J. Phys. B 44, 154008. 



19 



